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The phase diagram of the copper-oxide chain as a function of density and 
the nearest-neighbour Coulomb interaction, V, is determined. Phase separa- 
tion takes place above a critical value of V when the Cu 2+ — > Cu + valence 
fluctuations dominate. In the proximity of the phase separation boundary 
the superconducting correlations are the most divergent. We identify the 
parameter regions where the Luttinger Liquid theory applies and calculate 
the contours of the charge-charge correlation exponent K p . We show that 
anomalous flux quantization occurs as K p diverges. At n = 1, for V > t, a 
gap opens in the spectrum and the ground state is a charge-density wave. 
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The high temperature superconductors are difficult to understand quan- 
titatively because of the failure of mean field theories and perturbation meth- 
ods to give reliable results both in reduced dimensions and in the presence of 
strong electronic interactions. Under these circumstances exact calculations 
on finite size systems can provide valuable information. In particular, the re- 
sults of Luttinger liquid theory enable thermodynamic correlation functions 
to be deduced for one dimensional systems via exact calculations on finite 
size chains. In this letter we solve exactly finite length copper-oxide chains, 
up to a maximum of 12 sites, using the Lanczos procedure with the antici- 
pation that the copper-oxide plane will show some features in common with 
the copper-oxide chain [p]]. Our purpose is to study the phase separation 
and superconducting instabilities arising from the nearest neighbour copper- 
oxygen Coulomb repulsion. We identify the phase separation boundary of 
the extended Hubbard model via a Maxwell construction. Then we identify 
the region of possible superconductivity using the results of Luttinger Liquid 
theory and by the onset of anomalous flux quantization. 

Our model for the copper-oxide chain consists of two atoms per unit cell. 
Neglecting the oxygen-oxygen hybridization and considering the Coulomb 
interaction up to first-nearest neighbours, the copper-oxide chain is described 
by the two band model Hamiltonian, 
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where 

H t = -t J2 (dip j( r + h.c.). (1) 

<ij>a 

% and j are copper and oxygen sites respectively, < ij > represents nearest 
neighbours and the operator d\ a (p]- ff ) creates a Cu (O) hole with spin a. 
A is the charge-transfer energy, Ud (U p ) is the copper (oxygen) Coulomb 
repulsion, N s is the number of sites, N p is the number of particles, and V 
and t are the copper-oxide Coulomb repulsion and hybridization, respectively. 
The Hamiltonian has been defined so that it is invariant under the 'particle- 
hole transformation' d' ia — >■ d ia) Pj a — ► Pj a and A — > Ud — (A + U p ), i.e. 
E(JVp, A) = E(2iV s - JV P , U d -(A+ U p )). 



In this letter we consider parameters which are chosen to fit, as far as 
possible, with the photoemission data of high temperature superconductors. 
Thus Ud ~ 9t, A ~ 2 — 3t and U p small, where t ~ l.beV [§]. Hence the 
region 0.5 < n < 1.0 corresponds to Cu 2+ — > Cu + valence fluctuations, 
while (because of the 'particle-hole' symmetry of equation (1)) the region 
1.5 > n > 1.0 corresponds to Cu 2+ —* Cu 3+ valence fluctuations. As we are 
interested in the role played by the nearest neighbour interaction, V, this 
will be left as a free parameter. 

For most regions of the phase diagram the system is metallic, for which 
Luttinger Liquid theory can provide a valuable insight. However, for some 
parameter regions in the doping range 0.5 < n < 1.0 the system is unstable 
with respect to phase separation. This occurs at a finite value of the nearest 
neighbour Coulomb repulsion, and we identify this phase boundary before 
proceeding to a discussion of the metallic regime. 

Figure 1 shows the phase separation boundary as a function of V and 
hole density for U p — and U p = t. We set U& = 9t and A = 2t; the 
number of sites used was 12. The boundary for this region is found using 
the Maxwell construction. This is given by the points (ni,^) at which a 
tangent touches the curve of the ground state energy versus density. Other 
authors have used the point at which the compressibility, k, diverges (i.e. 
the point at which the ground state energy versus density becomes concave). 
These points, however, actually lie between n\ and n 2 and so could be a 
misleading indicator for when phase separation has actually occured. Figure 
1 also shows the line of k = oo. Using the Maxwell construction, the densities 
of the phase separated regions are then given by n\ and n 2 , and the fraction 
of the system occupied by n\ and ri2 at any given intermediate density n! 
is given by the relation mn\ + (1 — m)n2 = n', where m = f^ 2 ~" I . The 
tendency towards charge segregation can also be seen in the oxygen-oxygen 
charge correlation function ||, defined by 

C cdw (k) = -r 1 - £(n- < n >)(n i - < n >) e *( fl *-^), (2) 
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Where Hi = T,a4a C i,<r- 

The oxygen-oxygen charge correlation function (shown in figure 2) shows 
a peak at k — | (i.e. one over the system length) for both a density of 
n = 12 anc ^ n = if' This indicates a phase separation instability in the sys- 
tem. Phase separation occurs because, to avoid paying Coulomb repulsion, 
holes will tend to surround themselves with empty sites. This is because a 
repulsion, V, between a copper and a oxygen hole actually means an attrac- 
tion, — V, between a copper hole and a oxygen electron. The result is that as 
the doping is increased away from n = 0.5 holes will tend to move onto those 
oxygen sites which are surrounded by empty copper sites. However, where 
possible holes will stay on copper sites to avoid paying the charge transfer 



energy, A. The result is that the charge density becomes segregated and the 
Bloch symmetry breaks down [|J. 

The phase separated region is shifted as U p is increased into higher values 
of V, as can be seen in the inset of figure 1 . However, the phase separation 
boundaries do not depend strongly on U p provided that U p < Ud — A is sat- 
isfied. 

Having determined the region of the phase diagram for which a phase 
separation instability occurs we now turn our attention to the metallic phase. 
In the metallic region the long distance behaviour of the correlation functions 
can be determined via the aid of Luttinger Liquid theory. Before discussing 
our results we summarise some of the important aspects of Luttinger Liquid 
theory. 

Luttinger Liquids belong to the universality class of models described 
by conformal field theory with central charge, c = 1. This theory is useful 
because it gives relationships between the correlation exponents and a few 
simple properties of the finite system. Initially it was shown that models 
with only one degree of freedom, such as spinless fermions or bosons, had 
c = 1. Applying the same ideas to the Hubbard model proved difficult since 
conformal quantum field theory deals with Lorentz invariant systems only; 
i.e. it requires all gapless excitations to have the same velocity. For the 
Hubbard model this is true only at half filling, as then there is a gap for 
charge excitations. 

Recently it has been shown that conformal field theory can be applied to 
models with more than one gapless excitation if the excitations are decoupled 
in the low energy regime. In the low energy regime this is true of the Hubbard 
model f|. In such a case each degree of freedom has an associated operator 
algebra, which is again characterized by a central charge c = 1 . Provided that 
there are no phase changes as the interactions are increased it is reasonable to 
assume that Luttinger Liquid theory applies to lattice quantum models in the 
strong coupling regime. In fact, this assumption can be tested a posteriori 
by showing that the predictions of Luttinger Liquid theory are internally 
consistent. This is the assumption adopted in this paper for the extended 
Hubbard model. 

At long distances the charge-charge and spin-spin correlation functions 
are given by |J, 

. . , rt . K p cos (2kfx). o /2/ \ , cos (Akfx) 

{7TX) 2 X { - 1+K p> X 4iK p 

and 

< S z (x)S z (0) >= 7 1 ~ 2 + B 1 C °\ i ? k J* ) lnV\x), 

respectively The superconducting correlations are similar to the spin-spin 
correlator with K p being replaced by 1/K p . Hence, K p equals 1 for free 



fermions with spin. K p < 1 means that the charge and spin correlations 
dominate and K p > 1 implies that the superconducting correlations are the 
most divergent. 

Conformal field theory relates bulk quantities such as the compressibil- 
ity, k, and the Drude weight, D c , with the charge-charge correlation expo- 
nent. Where the system behaves as a Luttinger Liquid, the exponent K p 
is related to the compressibility by K p = ^nv c n 2 and to the Drude weight 
by K p = D c tt/v c (where v c is the charge velocity). When analytic solu- 
tions are not available these bulk quantities [0, || can be evaluated using 
exact diagonalization techniques . The discrete compressibility is given by 
N 2 k = 4/[e(N p + 2) + e(N p - 2) - 2e(N p )], where e(N p ) is the ground state 
energy per site. To calculate the Drude weight the chain is threaded with a 
dimensionless constant flux, <fi. This is done mathematically by the opera- 
tor transformation, c\ a — > c\ a e ^ . The finite-size equivalent of the Drude 
weight is then D c = [e(0o + 8(j>) — e(4>o)]/5(fi 2 , where 5(fi is an infinitesimal 
flux (~ 0.01) and 0o is a background flux which removes accidental degen- 
eracies of the ground state and guarantees a zero — current ground state. 
For example, O is (periodic boundary conditions) for N p = 4m + 2 and ■£- 
(anti-periodic boundary conditions) for N p = 4m with m an integer [10 . 



The charge velocity, as well as the spin velocity, can also be obtained 



from finite cluster calculations of the response functions [11]. In fact, they 



are related to the lowest frequency pole uji of the charge-charge and spin- 
spin dynamical correlation functions by vi = -^-, where q min = |p, N u 
is the number of unit cells and I refers to both charge and spin degrees of 
freedom |T3]. One can therefore directly determine both velocities without 



adiabatically following the relevant low-energy excitations, as has recently 
been done for the t — J model J7| . 

The phase diagram, figure 1, shows the contour of K p = 1, and the 
region for which K p exceeds 1. The expression used to calculate this contour 

was K p = TcnJ^j 1 . Notice that for low doping K p exceeds 1 within the 
phase separation region given by the Maxwell construction, and as doping is 
increased the curve K p = 1 moves out indicating that the superconducting 
correlations begin to dominate well within the conducting region. The points 
on the contour corresponding to the densities n = ^ and n = y§ are found 
from the 12 site chain, the point at n = yjj is found using a 10 site chain 
and the point at n — | is found using an 8 site chain. At high densities the 
point of Kp exceeding 1 is associated with the appearance of anomalous flux 
quantization. Anomalous flux quantization is the condition that the ground 
state energy is periodic in 0, with period tc, (i.e. the flux is quantised in 
units of tt). For the attractive Hubbard model, which is known to exhibit 
singlet superconductivity, the ground state energy satisfies 




£o(0 + 7r)-£ o (0) = ^(y) exp[-^), (3) 



where ( is a length associated with a gap of spin excitations and A depends 
on details of the model |T3| . As well as exhibiting anomalous flux quantiza- 
tion, the superconducting state is defined by the Drude weight (or superfluid 
density) remaining finite in the thermodynamic limit. For n = y| and yjj 
anomalous flux quantization appears almost exactly at the point that K p 
exceeds 1, the result being slightly better for n = ||. For lower densities this 
relationship is not so clear. 



We have also performed finite size scaling [J i] using a 6 and 12 site chain 
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at n = |. Using these two chains we are able to calculate the parameter £ of 



equation (3) and the ratio £ "\^ S= J for different values of V. We are then 
able to determine whether this is in the asymptotic region of equation (3) 
(i.e. the region of relatively large N s and small £). A small value of the ratio 
of the Drude weights in the asymptotic region will indicate that D c = in 
the thermodynamic limit, and that the true ground state is non-conducting. 
Anomalous flux quantization, without conductivity, is the signal for a pairing 
but non-superconducting state, such as a phase separated state or a charge 
density wave. When V — At (a point well within the phase separation region), 
C ~ 2 so this is in the asymptotic region, but the ratio g ,^"1 6 ~, ~ 0.1. Such 
a large reduction in Drude weight with such a small ( indicates an insulating 
state. Anomalous flux quantization first occurs when V ~ 2.7t for a 12 site 
chain at this doping. This indicates that at V = 4t there is a definite phase 
separation instability. 

For all other values of doping, namely n < 0.5 and n > 1, the effect of 
nearest neighbour repulsion is to suppress the superconducting correlations. 
For the densities n < 0.5 and n > 1.5, the phase diagram rather resembles 
the corresponding phase diagram of the t — U — V model ||. 

Finally, we consider the effect of the nearest neighbour repulsion at half 
filling, n — 1. At this filling the effect of increasing V is to force the system 
into an insulating phase. This can be seen from the graph of the Drude 
weight, the inset of figure 3, which vanishes at V ~ 1.6t. The driving force 
for this change is illustrated by the charge correlation function shown in fig- 
ure 3. It shows that as V is increased there is a sharp rise in the correlation 
function at k = ir (a wavelength of one unit cell). This corresponds to charge 
being pushed from copper sites onto oxygen sites, resulting in an insulating 
charge density wave. 

In conclusion, we have established the phase diagram of the copper-oxide 
chain using finite cluster calculations. In the region in which Cu 2+ — ► Cu + 
charge fluctuations dominate (i.e. A << U<i) we find that the coupling to 
the nearest neighbour repulsion leads to a region of phase separation for 
a sufficiently large interaction (typically V ~ A). In the proximity of the 
phase separation region and large enough density superconducting correla- 
tions are the most divergent. From our results we find that the curve K p = 1 
crosses the phase separation boundary at n ~ 0.72. At the densities n = ^ 



and j| the divergence of K p is accompanied by anomalous flux quantization. 
The contour K p = 1 has been calculated near the metal-insulator transition 
(n = 0.5) with only a small number of doped holes. The trend for the su- 
perconducting region to disappear asn-> 0.5 may be a consequence of this 
fact. Close to the phase separation boundary the superconducting correla- 
tions are precursory to phase separation. At weaker coupling the role played 
by charge transfer excitons [[jl| is still to be understood. For all other pa- 
rameter regions superconducting correlations are suppressed. At n — 1 the 
model shows a transition from a metallic state to a insulating charge density 
wave. 
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Figure 1: The phase diagram of the copper-oxide chain as a function of 
V and hole density for Ud = 9£, U p — and A = 2t, showing the phase 
separation boundary, the contour of K p —1 and the contour of infinite com- 
pressibility. There is also shown the onset of anomalous flux quantization for 
n = y§ and ^. The inset shows the phase separation boundary for U p — 
and t. 



Figure 2: The oxygen-oxygen charge correlation function for n = ^ and 
y|. This shows a peak at k — 7r/3 for V = 2t, indicating a phase separated 
state (U d = 9t, U p = and A = 2t) 



Figure 3: The Drude weight as a function of hole density for Ud = 9t, 
U p — and A = 2t, showing the opening of the gap at n — 1. The copper- 
oxygen charge-charge correlation function for a 12 site chain at n — 1 is 
shown for V = 0.5t and 2t. The density of electrons on the oxygen atoms 
rises from n = 1.265 to n — 1.797 in this range. 



